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VIRIAL EXPANSION BOUNDS THROUGH TREE PARTITION 

SCHEMES 

SANJAY RAMAWADH AND STEPHEN JAMES TATE 


Abstract. In this paper, we use tree partition schemes and an algebraic ex¬ 
pression for the virial coefficients in terms of the cluster coefficients in order to 
derive upper bounds on the virial coefficients and consequently lower bounds 
on the radius of convergence of the virial expansion T^vir- The bound on the 
radius of convergence in the case of the Penrose partition scheme is the same 
as that proposed by Groeneveld and improves the bound achieved by Lebowitz 
and Penrose. 


1. Introduction 

In this paper, we derive virial expansion bounds for stable potentials, based on 
the tree partition scheme of Penrose |23] . These bounds offer an improvement on 
the radius of convergence of the virial expansion T^vir when compared to those of 
Lebowitz and Penrose m and are actually the same as those proposed by Groen¬ 
eveld m Chapter IV, Section 3.7]. Towards achieving these bounds, we introduce 
a new approach to understanding the virial coefficients as a weighted sum over tree 
graphs. This offers an alternative approach to virial coefficients than those made 
in the literature. 

The virial expansion provides a description of the relationship between macro¬ 
scopic properties of a system of particles, specifically pressure and density. The 
virial expansion may be viewed as a generalisation of the ideal gas law VP = NkT. 
V denotes the volume in which the system is contained, P is the pressure of the 
system, N is the number of particles in the system, k is Boltzmann’s constant and 
T is the temperature of the system. 

The idea of refining the ideal gas law based on an expansion of pressure as a 
power series in density was proposed by Thiessen |31j and was first undertaken 
by Kamerlingh Onnes in 1901 m- An important reference for cluster and virial 
expansions is the book of Mayer m, who successfully derived the interpretation 
of the virial coefficients in terms of weighted two-connected graphs. Much work 
was done on the virial expansion in the 1950s and 60s, including: the develop¬ 
ment of Kirkwood-Salsburg Equations El; the graphical approximation approach 
of Groeneveld uni; deriving inequalities for the thermodynamic functions mm 
and bounds obtained from using Lagrange inversion on the cluster expansion m- 

In the literature, virial expansions have undergone a recent revival. There are 
many different approaches to the virial expansion. These include: a modification of 
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the Lebowitz-Penrose argument [2^; derivations through the Canonical Ensemble 
[24l [19] ; and using differential equations for functions related to the pressure with a 
uniform repulsive potential [3]. There is also the recent paper of Morals, Procacci 
and Scoppola m for the Lennard-Jones gas. Furthermore, the recent work of 
Jansen provides physical interpretations for the different radii of convergence, 
emphasising cases in which the virial expansion converges for a larger range of 
densities than the cluster expansion. 

A tree partition scheme for connected graphs was first provided by Penrose [23] . 
Sokal [2B] and Fernandez and Procacci [5] give a more general version and apply 
it to the Potts model and hardcore gases respectively. This paper develops this 
partition scheme for two-connected graphs. 

Providing a weighted tree expression for the virial coefficients should also have 
consequences for numerical results that have been obtained. The number of tree 
graphs grows at a smaller rate than those of two-connected graphs. Indeed, by 
Cayley’s formula, we have trees on n vertices compared to at least 2 ( 2 )“" 
two-connected graphs on n vertices. The lower bound arises from considering the 
graph corresponding to a cycle on n vertices, which is certainly two-connected. We 
may add or not any of the remaining ( 2 ) — n edges in the graph to give some 
subset of two-connected graphs. The work of Ree-Hoover [25] established a method 
of reducing the number of two-connected graphs that one needs to consider by 
introducing a factor 1 = hij — fij for every missing edge in each two-connected 
graph and finding the cancellations that arise from this. This has led to Monte- 
Carlo simulations to give many virial coefficients in the work of Clisby and McCoy 
[5]. A practical application of these expressions would be to evaluate the tree 
versions provided by the expression of virial coefficients in terms of weighted trees. 

In future work, we hope to apply this method to other partition schemes in 
particular to the scheme indicated by the notion of internally and externally active 
edges of a matroid [28l [2], which is given explicitly in m- Furthermore, the 
hardcore cases in [9] will give good further examples to apply the tree partition 
scheme to. 

The paper is structured so that in section [2] we give a derivation for the re¬ 
lationship between cluster and virial coefficients emphasising their combinatorial 
properties. The results of this paper are summarised in section [3] and we prove that 
this result is the same as that proposed by Groeneveld in section[5]. In section[6l we 
provide the general framework in which we understand the tree partition schemes 
and the important properties. We apply this to the Penrose partition scheme in 
section [7] for positive and stable potentials and present the conclusions and further 
work in section [5] 


2. The Classical Gas Model 


The results for this paper are for the classical many particle system. 

The canonical partition function for a fixed number of particles is given by: 



( 2 . 1 ) 


where Un is the interaction potential. We make important assumptions on the 
interaction potential in order to achieve the results, which are given in subsection 
\TT\ In subsection 12.21 we indicate the connection between cluster coefficients and 
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weighted connected graphs and in subsection 12.31 we derive the algebraic relation¬ 
ship. 

2.1. Assumptions on the potential. 

Assumption 1 (Potential). The N-particle interaction potential: Un{xi, ■ ■ ■ ,xn) 
may he written as the sum of pair-potentials: 

Un{xi,-■ ■ ,xn) = ^ ^{xi,Xj). ( 2 . 2 ) 

Furthermore, we assume that the pair potentials ^{xi,Xj) are central, that is, they 
only depend on the distance from Xi to Xj. 

Definition 1 (Stability). A pair potential is called stable if for all N and all 
collections of locations (cci, • • • ,Xn) S , we have the following inequality: 


^ ^{x„Xj)>-BN. (2.3) 


Definition 2 [C{(3) - ‘temperedness’). The ‘temperedness’ function C{j3) plays an 
important role in the estimates of convergence through tree partition schemes. 



g-/3$(0,a:) 



(2.4) 


If the expression (EH) is finite then the potential <i> is called ‘tempered’, which we 
assume for bounds involving C{fi). 


2.2. The Classical Gas and Weighted Generating Functions. We use the 

finite-volume grand canonical partition function for a classical gas: 

OO 

2a(z) = ^ (2.5) 

N=0 

where, given the assumptions on the potential, ^n.a may be expressed as: 

We have the relationship /3 Pa = j^lnSA, which gives us an expansion for the 
finite-volume pressure. 

The cluster and virial expansions for the finite-volume pressure are expressed as: 


,3P., = f; 

n\ 

n—l 

(2.7) 

= f; 

n\ 

n—l 

(2.8) 


The variable z = is the activity or fugacity of the system and gives control 
over the average particle number in the grand canonical ensemble. The variable 
fi is the chemical potential and fi = -^ is the inverse temperature, p denotes the 
density of the system of particles. The thermodynamic pressure p is then found by 
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taking the limit A t in the sense of Fisher. The cluster and virial 

the thermodynamic pressure are expressed as: 

expansion for 

P= 

n\ 

n—l 

(2.9) 


(2.10) 

The fugacity and the density are related to each other as follows: 


dp bn n 

dz ^ (n — 1)! 

n—l ^ ' 

(2.11) 


Derivations and expressions of the coefficients 5„(A), /3„(A), bn, and /3„ in terms 
of connected and two-connected graphs are found in the book of Mayer [181 Chapter 
13], 

Mayer’s trick is to introduce fij := exp{—/3^{xi,Xj)) — 1 and to rewrite the 
integrand as: 

exp(-/3 exp(-/3$(a;,,a;j)) 

= n (i+/*j) 

= X! n 

SGS[N] eGE(g) 

where t/[-/V] denotes the collection of all graphs on the vertex set [N], 

We may define the graph weight: 

w{9)-= n (2.12) 

e£E(g) 

then the canonical partition function may be written as: 

2’A,A = ^n(/ '^^9) (2-13) 


and we realise that the grand canonical partition function may be written as the 
weighted generating function of graphs. One may also use the weights including 
the integrals as: 



(2.14) 


Mayer derived that the pressure can be written as the weighted exponential 
generating function of connected graphs, that is, for the weight: 


AfR"^ |A| 


(2.15) 


We have that: 

bn= ^( 5 )’ 

geC[n] 

where C[n] denotes the collection of connected graphs on n labels. 


(2.16) 
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2.3. Deriving the Algebraic Relationship. It is possible to relate the virial co¬ 
efficients /3n to the well understood cluster coefficients bn- This is done by inverting 
(12.111) . so as to obtain an expression for the fugacity as a power series in terms of 
the density, and then substituting this power series into (12. 9|) . Lagrange inversion 
[6] gives the following: 


Pn+l = n\[z"] 



(2.17) 


Using this expression, we can find a nice combinatorial expression of the /3„ in 
terms of the bn- 

Define, for any a; G R and any nonnegative integer k, the falling factorial as: 


fe-i 

(x)- := x{x — 1) • • • (x — fc -I- 1) = {x — i)- (2-18) 

i=0 

If X is a positive integer, we may also write (x)- = and = (J;). We 

use this last equality to extend the definition of binomial coefficients. For x G K 
and any nonnegative integer k, we define 



(2.19) 


Let Xi, X 2 ,... be an infinite number of variables, then the partial Bell polynomials 
Bn,k {.xi,X2, - - - ,Xn) are dehned by: 


exp I u ^ 2:™^ I = 1 + X! (a^i,... ,x„) j . (2.20) 


m—1 


n—1 


n\ 


\k^l 


There exists an explicit expression for the partial Bell polynomials: 


^n,k (^1; ^2; • ■ • ; ^n) 


E 

(ki . ,A:n):fci >0 
Ei fci=fc^Ei iki=n 


n\ 





( 2 . 21 ) 


For more known results on Bell polynomials we refer the reader to [6] . Lastly, we 
define the potential polynomials. Let xq = 1, let xi,X 2 ,... be an infinite sequence 
of numbers and let r be any number. The potential polynomials (xi,..., x„) 
are defined by: 


E"^-^ =l + ^pW(xi,...,x„)-. 


( 2 . 22 ) 


\n—0 


There is also the following explicit expression for the potential polynomials: 


(xi,..., x„) = y^^{r)-Bn,k (xi,..., x„). (2.23) 

fc=i 


We return to (I2.17|) . Since ^ = E^i (Ei)! ^" ^ 
it follows from (12.221) and (12.231) that: 
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^n +1 = n\ - j - 

n\ 

n 

= y^(-n)-Pn,fc (62, ■ • ■ , bn+l) 
k=l 

= (~k) (^ 2 ,..., 6 .+ 1 ). (2.24) 

fc=i ^ 

3. The Main Result 

The main result of this paper is the improved radius of convergence of the virial 
expansion. En route to this improved radius of convergence, we obtain a new 
expression for the virial coefficients in terms of weighted tree graphs. This new 
expression is given in section [71 

For a stable interaction, we let u := exp(2/3P), where B is the stability constant. 
We have the following bound on the radius of convergence of the virial expansion, 
Pvir; achieved through using tree partition schemes for connected graphs, arising 
from the Penrose partition scheme. The proof for the following theorem is given in 
section [71 The function Ti(z) is the following generating function: 

°° ^n+l 

n—1 ^ ' 

Theorem 17.21 For a classical gas with a stable potential, the virial coefficients 
have the following upper bound: 

^ / l + lTi(uc(u)) 

(n + 1 )! “ n \ c{u) 

where c{u) is the smallest positive root of 

c{u)T[{uc{u)) - Ti{uc{u)) = 1. 

u 

This result is precisely the same as the one proposed by Groeneveld m Chapter 
IV, Section 3.7], which is proved in section |5l 

4. Comparison of the bound to previous bounds and application to 

PARTICULAR MODELS 

We have the exact equation (ED for the virial coefficients, which may be used 
to calculate the coefficients in particular examples. One just has to calculate the 
weights of non-splittable Penrose trees. 

For stable potentials we have a bound that is an improvement of Lebowitz 
Penrose[T5] that we can see in in figure [T] In particular for positive potentials, 
we obtain 7?.vir > 0.237961C'(/3)“^, whereas the Lebowitz and Penrose bound gives 
7^vir > 0.144766998C(^)-b 


(3.2) 

(3.3) 
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Figure 1. A comparison of the lower bound for the radius of 
convergence of the virial expansion obtained by Lebowitz-Penrose 
with that given by Groeneveld 

5. Penrose is the same as Groeneveld 

In m, Groeneveld obtains, for stable potentials, the bound TZvir > a{u)C{P) 
where a{u) is the smallest positive root of method relies 

upon finding recursion relations for graphs, expressing these as a system of differen¬ 
tial equations in three variables for multivariable generating functions and solving 
these equations in order to obtain an expression for the virial coefficients. Here we 
indicate how our result actually gives the same. 

In order to show that the two bounds are equivalent, we use the fact that the 
rooted tree generating function is the inverse of s i-A se“® , which can be found in 
[7]. The function Ti{z) as defined in (I7.12|) is shown to be equal to Tpen,i{z) in 
section[3and we use in this section two identities proved in sectionOfor Tpen,i(-7), 
but we write these in terms of Ti{z) to maintain consistency with the notation of 
the previous section. 

Lemma 5.1. The derivative of the generating function Ti{z), evaluated at se~'^ 
satisfies the following identity: 

T[{se-^) = (5.1) 

1 — s 

Proof. We use the relation (I7.30|) : T{{z) = ■ If we change variables z M- se~^ 

and use that T*(se“®) = s, then we easily obtain: 

Tiise-^) = (5.2) 

I — s 

□ 


Lemma 5.2. The generating function 7i(z) evaluated at se ® satisfies: 

Ti{se-^) = I -e-". 


(5.3) 
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Proof. We know from (j7.28ll that Ti(z) = 1 — making the substitution z i-A 

se“®, we immediately get: 

ri(se-") = 1-e-b (5.4) 

□ 

Proposition 5.3. The lower bound of the radius of convergence of the virial ex¬ 
pansion: 


Tevir > Cil3) 


-1 


z{u) 


1 + ^Ti{uc{u)) 


(5.5) 


where c(u) is the smallest positive solution to uc{u)T[(uc(u)) —Ti{c(u)) = 1 is equal 
to a{u)C{l3)~^, where a{u) is the smallest positive solution to a(u)e 


— a{u) _ 


(l+ii)e ' 

Proof. Reparametrise uc(u) = t(u)e“*^“^ and substitute into the defining equation 
for c{u)\ 

t(u)e-‘(“)T((t(u)e-*(“)) - ri(t(u)e-*(“)) = u. (5.6) 

We use dSI]) and ()5.3p to simplify the expression to: 


1 — t{u) 


- 1 + e 


-t(ti) _ 


-t(«) 


1 — t{u) 

We transform a{u) = 1 — t{u) and obtain: 

(1 + u)a{u)e = e“^“\ 

1 


= 1 + M. 


(5.7) 


a{u)e 


— a(it) _ 


(1 + u)e' 

For the bound, we make the same substitution: 


"^Vir > 


and use (15.3|) to obtain 

We use dEZl) to obtain 


T^Vir > 


u + ri(t(ii)e“*(“))' 
t{u) 


(1 + u)e*(“) - 1' 


'^Vir > 


t{u) 


1 _ 1 
1 —i(ii) 


= 1 — t(u) = a(u) 


as required. 


(5.8) 

(5.9) 

(5.10) 

□ 


6. Partition Schemes 

The improved bounds of this paper rely upon partition schemes for connected 
graphs in terms of tree graphs. A helpful reference for partition schemes and their 
application in the cluster expansion of hard-core models is found in [5]. 

We define [n] := {1, • • • ,n}. The collection of all connected graphs, respectively 
tree graphs, on the label set [n] is denoted by C[n], respectively T[n]. 

We define a partial order on C[n] by bond inclusion: G <G E{G) C E{G). 
For G < i7, we define the set [G,H] = {K\G < K < PI}. We consider partitions 
of C[n], indexed by t G T[n], comprising of blocks of the form [r, i?(T)], where 
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R : T[n\ C[n]. The mapping R is called a partition scheme. More formally, we 
have the following definition: 

Definition 3 (Partition Scheme). A partition scheme for a family C[n] of connected 
graphs is any map R : T[n] C[n\, t ^ R{t), such that: 

i) E{R{t)) D E{t) for all r G T[n], and 

ii) C[n] is the disjoint union of the sets [t, R{t)] for r G T[n\. 

The advantage of a Boolean partition scheme is that it allows us to rewrite the 
connected graph sum as: 


E n /«= E n n (^.i) 

GeC[n]eeE(G) reCfra] eeEfo) eeE(R(T))\E(T) 

This section outlines the general method of achieving bounds from partition 
schemes. In the next section, we apply this to the Penrose partition scheme. 

This approach requires, firstly, a well defined mapping from products of trees 
arising from cluster coefficients in (12.241) to trees on the vertex set \n + 1] such that 
the weights are preserved. This is achieved through the product structure. It is 
then necessary to understand the combinatorial factors that arise from this many- 
to-one mapping. After achieving an exact expression, we make appropriate bounds 
for the tree weights, depending on whether we assume positive or stable potentials. 

6.1. Faithful Mergings of Trees. In expression (12.241) we have an expression for 
the virial coefficient as a many variable polynomial in the cluster coefficients. This 
polynomial is explicitly given by: 


Pu+l = ^ 


k=l 


—n 


kl 


{kl 

Eifc 


E 

,...,fcn):fci>0 

i = fciEi ik^=Ti 


nr= 


n 


^ 1—1 


H-\-l 


( 6 . 2 ) 


Equation (IQ) shows how we can relate the cluster coefficients to a sum over 
trees using any partition scheme. Now, consider any monomial ]([j appearing 
in (IQ) . We want to find an analogous expression for this monomial in terms of 
trees on n + 1 vertices. 

Using any partition scheme, the monomial can be related to a sum of 

products over trees. Any term in this latter sum is a product which involves a 
collection of k trees, say {ti, ... ,Tfe}. We want to associate to this collection of k 
separate trees a single tree on n + 1 vertices. This is done by merging trees. 


Definition 4 (Merging of trees). Let I be a set of labels and let {ti, ... ,Tk} be 
a collection of I-labelled trees. Let G be the graph with vertex set L and edge 
set G Ti for some 1}. The edges come with multiplicity rn{{i,j}) = 

\{1\E{ti) ^ Then G is called the merging o/{ri,..., t^}. 


It is important to note that the merging of a set of labelled trees heavily depends 
on the labeling of the trees. In general, a merging may contain multiple edges 
between a pair of vertices, loops and cycles. It may even not be connected at all! 
See Figure I for an example. 

If a merging of a set of labelled trees is also a tree, we say that the set of labelled 
trees is properly labelled. 

Denote by Vi the label set for the tree Ti arising from a merging. 
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Figure 2. An example of ill behaviour of mergings. On the left- 
hand side we have three trees labeled in {1,2,3,4, 5}, and on the 
right-hand side we have the merging of these trees. 


Proposition 6.1. A necessary condition for the merqinq to he proper is that \Vi O 
V,\<lforain^je[k]. 

Proof. If |Vi n FjI > 2, this means we have (at least) two vertices say I and m 
in common between Ti and Tj. The first problem is whether we have a repeated 
edge in common between the two graphs, if we do then it is not properly labelled. 
Otherwise we have two disjoint paths from I to m arising from the path in Ti and 
the disjoint path in tj since they are trees. This means that we will have a cycle 
in the merged graph and it therefore won’t be a tree. Hence |Vi n Vjj < 1 is a 
necessary condition. □ 

We now assume that |Fi n V, | < 1 for all i ^ j G [k] and provide a necessary and 
sufficient condition for the merging to be proper. 

Definition 5 (Merging Graph). For a qiven merging with label sets (l{)ie[fc]: 
define the merging graph A4((Vi)ig[fe]) through its edge and vertex sets: 

F(A4((F0^6W)) := [fc] U ( (J (F. O F) x {*} 

F(A4((F).6[fc])) := {{^, (F n F,*)} {j, (F n F,*)}| n F ^ 0}, 

where 0 x {*} := 0. 

Proposition 6.2. Given that the vertex sets F satisfy |F O F| < 1 for all i ^ j G 
[k], a merging o/ (ti, • • • ,Tfc) is properly labelled if and only if the merging graph 
A4((F)*e[fc]) is a tree. 

Proof. To prove the proposition we use the two defining properties of a tree: that 
it is connected and acyclic. 

The merging graph A4((F)ie[fc]) is connected if and only if the merged graph G 
is connected since there is a path between I G Ti and m G Tj in graph G if and only 
if i and j are in the same connected component of At((F)iG[fe])- 

We now show that there is a cycle in the merging graph Af((F)ie[fc]) if 
only if there is a cycle in the merged graph G. If we have a cycle in At((F)ig[fc]), 
then replacing each vertex i in the cycle by the corresponding path in Ti between 
its two neighbours in the cycle, we end up with a cycle in G. The cycle in G 
can be transformed into a cycle in AI((F)ie[fe]) by reducing each tree path to the 
corresponding vertex as no cycles will appear in an individual tree. □ 
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Remark 1. We realise here that the graph Ai{{Vi)i^[k]) is given a similar struc¬ 
ture to the block outpoint tree used to understand the two-connected components of 
connected graphs. The vertices at which we attach two trees are like the outpoints 
and each individual tree is treated like a block. 

Now, let R be any partition scheme and consider any product Hi ^^+i appearing 
in As mentioned earlier, this product can be written as a sum of products 

of k trees. Consider one such product and say it involves the collection of trees 
{ti, ..., Tfc}, then we can associate a tree on n + 1 vertices, say t, to this collection 
by properly labeling the trees ri, ..., t^,, where the labels are chosen from [n + 1]. 
The tree r is then simply the merging of {ri ,... ,Tk}. Since the k smaller trees 
have = n + fc vertices in total, it follows that we have to identify fc — 1 

vertices. While any proper labeling with fc — 1 identifications will do, we require 
that the merging respects the partition scheme R in the following sense: 

Definition 6 (Faithful merging). Let R be any partition scheme and let {ti, ..., T/t} 
be a collection of properly [n-\-l]-labelled trees. An R-faithful merging o/{ti, ..., Tk} 
is a labelled tree r such that: 

i) E{t) = ujLiA(ri), and 

ii) E {R{t)) \ E{t) = Uti {E {Rin)) \ E{n)). 

We call {ti, ..., T}~} a splitting of t. 

We define the weight of a tree a;(r) by: 

w(r) := n n (l + A)- 

eeE(T) eGE(R(r))\E(r) 

The importance of a faithful merging follows from the following 
property. 

Lemma 6.3 (Product Structure of Faithful Mergings). Let {ri,...,rfe} be given 
and let t be a faithful merging of this given set of trees. We have the identity: 

k 

]^w(tO = a;(r). (6.4) 

Proof. From property i) of a faithful merging, we have that: 

n n /•= If /. 

i—'i-eGE{Ti) eGE{T) 

Property ii) of the faithful merging gives that: 

k 

n n ( 1 +/^)= n 

i=l eeE{R{Ti))\E{Ti) eGE(R(r))\E(r) 

Taking the product of the two identities, it follows that: 

k 

If If If (1+A)= U fe n (1 + /-) 

i=leGE(n) eGE{R{Ti))\E{Ti) eeE{T) e&E(R(T))\E(T) 


(6.5) 


( 6 . 6 ) 


(6.7) 

□ 


(6.3) 

factorisation 
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6.2. Factorisation of Integrals. We wish to show that the full weight of a tree 
also factorises as we require. 

The full weight of a tree is given by 



( 6 . 8 ) 


Since we are assuming that the pair potentials are centred, the integrand de¬ 
pends only upon the distances between the particles. The thermodynamic limit is 
independent of boundary conditions and so we consider periodic boundary condi¬ 
tions. With periodic boundary conditions, we may make the change of variables 
from the spatial points to the differences. We may choose the change of variables 
to be {xi, ■ ■ ■ ,xn)'-^ ixi,X 2 — xi, - ■ ■ ,xn — Xi) and in this case all differences can 
be expressed in terms of the latter — 1 variables and so is independent of xi 
integral. The key point is the periodic boundary conditions, which allows for the 
change in coordinates to preserve the area over which we integrate. 

The factorisation of the integrals occurs analogously to the block factorisation 
of connected graph weights as explained in m- The graph acts in 

the same way as the block cutpoint tree, where the vertices corresponding to trees 
bear the analogy with blocks and the vertices corresponding to the common label 
bear the analogy with outpoints. First we give the definitions required to rigorously 
convey how to treat the factorisation of these integrals. 

Definition 7 (The Centre of a Tree). The eccentricity of a vertex v in a tree t is: 
e{v) = max{(i(u, w)|r(; G V{t)}. 

The radius of a tree is: radius^r) = min{e(u)|u G F(r)}. 

The centre of a tree is: {v G V{T)\e{v) = radius^r)}. 

Remark 2. A priori the centre of a tree may be a set of cardinality strictly greater 
than 1. We are working with a special class of trees for A4{{Vi)i^[k])> which are 
bipartite and importantly have all their leaves in a single set. This is sufficient to 
get a unique centre. 

The idea is that the vertex w for which the eccentricity of a vertex is attained is 
necessarily a leaf. If not then we always have a neighbour away from the original 
vertex which allows us to increase eccentricity by 1. So any intersection point has 
odd eccentricity and any tree vertex has even eccentricity and so the centre can 
contain only vertices from one of the two sets in the bipartite graph. The centre 
cannot contain two points of distance at least two from each other and hence can 
only be a singleton, so the centre of AI((I^)ig[fe]) is well defined as a single point. 

When we want to understand the factorisation, we want a well-defined way of 
deciding what label to take as the redundant integral in each of the tree factors. 
When we see how to do this systematically, we achieve the factorisation of the tree 
weights for the tilde generating function as well. 

In order to define the vertices over which we make our integrals independent, we 
dehne a digraph from the tree by orienting all edges away from the centre of the tree. 
In this case every vertex except the centre has a unique vertex pointing towards 
it. For each tree vertex outside of the centre, the intersection point neighbour 
that is the other endpoint of the edge pointing towards it is the vertex, which the 
integrand is taken independent of. Each intersection point outside of the centre 
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has also precisely one incoming edge which marks the tree for which this point is 
integrated over. 

We now consider what happens at the centre. If the centre is an intersection 
point, then this is understood as the overall point which is not integrated over. If 
the centre is a tree point, then we may choose any of its neighbouring intersection 
points to be the independent point and we are left with this single independent 
point. 

Hence we have that the tilde weights factorise. 


6.3. Establishing the Combinatorial Relationship. By combining equations 
dH]), (ESI), and EH) we find: 






Si ^i~^>Si 


_iL_n 






n\ 


h-\-i 


n 


ki 


Ei ki=k,^i iki=r 


nr=ifc.!(*! 


ifcf n^H-i 
2=1 




= E(T)^' E 


n\ 


>0 Ili=l'^*-Ud 


Ei iki = 


(ri,...,T|c) i —1 


where the sum is over Ti, which are spanning trees for the corresponding bi in the 
product on the previous line. The next point is to write this product of the weights 
of tree graphs as the weight of an i?-faithful merging. This involves assigning labels 
to the individual trees as described above. We then express it as: 


Pn+l — 



E 

(fej ,... ,kn )'ki'>0 
Y,i fci=fe)Ei ik^=n 


ir:=i 


E 

(tI,-" ,Tk) 


1 

N{ti, ■■■ ,Tfc) 


E 

TGT[n+l] 


(6.9) 

where in this case r is some R-faithful merging of (n, • • • , Tk) and fV(ri, • • • , Tk) 
is the number of such i?-faithful mergings. 

We want to express the latter sum as a sum on 7'[n + 1] and drop the condition 
that the tree graph is an i?-faithful merging of a given set of k smaller tree graphs. 
Note that this can be done, since any t £ T[n + 1] can be made its own i?-faithful 
merging by choosing a suitable labeling on its vertices. However, if a tree graph 
r has a splitting {ti, ..., t^}, then it also has splittings {fi, ... ,fj} with j < k 
and with the fi being mergings of distinct Tj, meaning that some trees will appear 
multiple times in the last sum in the right-hand side. 

We can overcome this problem by rewriting this sum as follows. We call a tree 
graph r Z-splittable if r has a splitting into I parts, but not into I + 1 parts. The 
number of Z-splittable tree graphs is a combinatorial factor which depends on the 
partition scheme being used. Combining this factor with the other combinatorial 
factors yields: 


I3n+i = ^C{n+l,l) ^ w(t). 

; = 1 T€Tln.+ l]: 

T is J-splittable 


( 6 . 10 ) 
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where 

i 

C{n + lJ) = J2 




—n 

k 


k\ V _ - _ V _i_ 

n” 1 


Ei ki=k,J2i iki=n 

with the sum over (ri, • • • , r^) being over splittings of r. 


( 6 . 11 ) 


Theorem 6.4. The virial coefficients may be expressed in terms of the exponential 
generating function of non-separable tree graphs according to a partition scheme R 
as follows: 


/3n+l 

(n + 1)! 


1 


n + 1 


^qR,iC{n + 1,1) 


1^1 


{Tn,i,Uz)y 

l\ 


( 6 . 12 ) 


where guy is a combinatorial factor denoting the number of ways to attach non- 
splittable graphs to make an l-splittable graph. 


This theorem is a rewriting of (16.101) . where we realise that an Z-splittable graph is 
uniquely defined by its I non splittable components, which arise from the generating 
function Tr, 1 ,^( 2 :) and the combinatorial factor ensures each Z-splittable graph 
is used only once. 


6.4. Obtaining Bounds from the Tree Partition Scheme. The final stage is 
then to make bounds according to the assumptions made on the potential. The 
method of obtaining such a bound is to give a uniform upper bound for the product 
of 1 + /-factors and then to calculate the integral of the |/|-factors per each edge. 
The value of the integral is C{f3) as defined in section [2l 
For stable potentials we have the following theorem. 

Theorem 6.5. For u := exp(2/3il), we have the following bounds for the virial 
coefficients, when we have a stable potential. 


Wn+i\ < C'(/3)”u"-i ^ \C(n + 1, mraAn + 1]| 


Z =1 


n -|- 1 ^' 


(u -I- 1)! 


1=1 


n 


(6.13) 


(6.14) 


This is a direct consequence of the following lemma and (I6.10L respectively 
theorem 16.41 for (16.131) . respectively (16.141) . 

Lemma 6.6. There exists a permutation of the vertex labels 11 such that when 
applied to the labels in the product: 

n (1+/-)’ ( 6 . 15 ) 

eeEiR{T))\E{T) 

we have the upper bound for r a tree on n + 1 vertices. 


Proof. We construct the permutation 11 required above, by using an important 
observation found in the book of Ruelle m about bounding 1 -|- /-factors. 

Recall that stability requires that for any subset S' C [n -|- 1], we have that: 

^{xi,Xj) >-/3B\S\. (6.16) 

{z.j}GS(2) 
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This implies that 

n il + kj)<exp{PB\S\). (6.17) 

{ijleSG) 

We realise that the square of (16.151) may written as: 

nn (1 + /*.,) <exp(2/35|5|), (6.18) 

ies jGS 

and hence that one of the factors in this product is less than exp(2/3i?). We call 
this label 15 . Therefore we have: 

Y[{l + f.s,j)<eMWB). (6.19) 

jes 

In order to bound the product 

(1 + /e); 

eeE(R(r))\Eir) 

we consider the vertex labelled n + 1 and its neighbours in the graph Gn+i, defined 
through its edge set E{Gn+i) ■= E{R{t))\E{t) and call this neighbourhood N{n + 
1). We let Sn+i = iV(n + 1) U {n + 1} and then relabel n + 1 i-A and 

I'Sn+i n + 1. We call this transposition ct„+i. 

This gives us that: 

n (1 + /-n+i(n+i).,) < exp(2/3il). (6.20) 

jGSn+l 

j^cr„+i{n+l) 

The permutation 11 is defined through a sequence of transpositions, which are 
defined inductively as below. 

For each stage in the induction, we retain the original labelling of the vertices in 
order to find which transposition is required to apply to obtain the bound. Given 
the relabellings for [i + 1 ,n + 1 ], we understand how to relabel i. In the graph Gi 
defined by edge set E{Gi) = E{R{t)) \ E{t) fl and vertex set [i], we consider 
the neighbourhood of i and call this N{i). We need to apply the sequence of 
transpositions generated so far to this set of ’original’ labels in order to understand 
how they are labelled at the current stage. We use the notation ttS := { 7 r(j)| j G S'}, 
where tt is a permutation and S is a set. 

We define 

( II ] ({0 U iV(i)) = (Ti+i ■ • •cr„+i({i} U A^(i)). 

\j=i+i J 

We then define the transposition 

o'i = n • 

\j=i+i J 

The required final permutation is 

n+1 

n — (Ji — (J3 • • • 

2=3 
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which is what we apply to our labels. 

This process provides the identity: 

n+1 

= n n (6-21) 

j=3 j&Si 

We apply this permutation 11 to the integrand since the final result is independent 
of the labelling of the particles. 

Each factor in this product is bounded from above by u := exp(2/3i3) and we 
have n — 1 of these. We thus have a factor in the virial expansion bound. □ 

In the following section we illustrate this method for the Penrose partition 
scheme. 


nl n (1 

^sGE{R{t))\R{t) 


fe) 


7. The Penrose Partition Scheme 


In 1967, Penrose [53] introduced a partition scheme in the sense of Section 131 
The scheme is as follows: 

Suppose we have a graph G with labels in [n + 1] and let r G T[n + 1]. For any 
vertex label i of r, let d{i) be the tree distance of the vertex labelled * to I and let 
i' be the predecessor of i i.e. d{i') = d{i) — I and {*',*} G E{t). We associate to 
T, the graph i?Pen(T) found by adding (only once) to r all edges {i,j} & E{t) such 
that either: 

PI d{i) = d{j) (edges between vertices at same generation,) 

P2 d{j) = d{i) — 1 and i' < j (edges between vertices one generation away.) 

The first result achieved is an expression for the virial coefficients in terms of the 
weighted trees. 


Theorem 7.1. The virial coefficients may be expressed in terms of weighted trees, 
with respect to the Penrose partition as follows. 

where 


Tpe„T.i(z) = 


.l\ 

n—1 Te7pen,l[’T.+ l] 


(7.2) 


As a consequence of the above writing of virial coefficients in terms of weighted 
trees, we obtain the following bounds for the virial coefficients and the radius of 
convergence of the virial expansion. 


Theorem 7.2. For a stable potential, the virial coefficients are bounded as: 
\l3n+l\ ^ C{I5Y A + Apen.l (mc('u)) 


< 


(n + I)! n + I 


z(u) 


(7.3) 


where c{u) is the smallest positive solution to 

uc{u)Tf^^^^{uc{u)) - Tpen,i{uc{u)) = u. (7.4) 

This provides the bound on the radius of convergence of the virial expansion as: 

c(m) 


T^Vir > 


1 + iTpen,l(nc(u)) 


cm- 


(7.5) 
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1 



Figure 3. An example of a Penrose-faithful merging. The labels 
in the vertices on the left-hand side are the original labels, and 
the labels by the vertices are the newly assigned labels. On the 
right-hand side we have the merging of these trees. 



7.1. The Combinatorial Factors. 


Lemma 7.3. The number of different Penrose-faithful mergings for a collection of 
trees (n, • • • ,Tk) is: 


N{ti, • ■ • ,Tfc) = fc! 


n\ 


(7.6) 


Proof. In order to assign labels to the individual trees in a monomial to obtain 
a Penrose-faithful merging, we first order the individual factors. This contributes 
kl to the N{ti, ■ ■ ■ ,Tk). Once the trees are ordered we assign new labels to these 
trees. The trees all come with their original labels from the expression for connected 
graphs. We respect the ordering of the labels. We assign the label 1 to the first 
tree, this is given to the vertex originally labelled 1. We then assign i labels to each 
tree with i \ vertices. This is done in }^\[i\)K ways. We assign these labels 
in their natural order to the vertices originally labelled in [2,i The vertex 1 

achieves a label from the preceding tree. From the preceding tree, we choose the 
vertex that is firstly the greatest distance from 1 in the orignal tree and if there are 
more satisfying this property, we choose the smallest. 

□ 


We realise that for an m-splittable tree we have ways of producing this 

tree from a product of k terms. We have m — 1 points at which to split the tree 
and we split at fc — 1 of these. The combinatorial factor in (I6.10|) is therefore: 



Proposition 7.4. We have the combinatorial identity: 



(7.7) 


(7.8) 
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Proof. If we let s = fc — 1 in dm), then we can recast the sum as: 


m —1 


E 



(7.9) 


We realise this is the ^ term of the generating function formed by multiplying 
(1 + and i ^ {i+zY ~ ’ '^hich means we need to calculate: 


1 


(1 + ^) 


n—m +1 


-(1 + ^)^ 




(7.10) 

□ 


Hence, C{n + 1, m) = (~l)"*(m) and we have the expression: 

A+i = i:(-ir(") i: iiir)- (7.11) 

771=1 ^ ^ reTpen.mh+l] 

Let rpen,i( 2 :) be the shifted exponential generating function for Penrose graphs 
that cannot be split. It is defined by: 


00 

T^enAz) := ^ IrPenaWlT^-^. (7.12) 

^2 - 1 )' 

The exponential generating function is shifted to make it easier to find the num¬ 
bers |7pen,mHI from the above series. 

To form an m-splittable tree, we need to attach m non splittable trees together 
in the way described above. We have that 

|rpe„.^Ni= ^ ^^z_!^]^|rpe„,i[* + i]|'= 7 , (7.13) 

where the sum above is over sequences (fci,-- - , fc„) satisfying ^ and 

iki = n. We can thus see this as the coefficient of in the product Tpen,i(z)'" 
and hence we obtain: 

irPen.mNI = (n - l)![z"-i](Tpe„.i(z))"*. (7.14) 

We note that this identity is true for the weighted generating function T-pen,i,Cb{z) 
as well due to the factorisation of the weights and hence 

Tpen.mA^) = (Tpen.l.ii ( 2 ) )™ • (7.15) 

7.2. The Generating Function rpen,i(2)- The main result of this section is the 
number of non splittable trees with respect to the Penrose partition. 

Theorem 7.5. The number of non-splittable trees for the Penrose partition is: 

|rPenT[n+l]| = (n-l)"-\ (7.16) 

where we understand as 1 to eover the ease n = 1. 

In order to prove this, we first require modifications of the following two theo¬ 


rems. 
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Theorem 7.6 (The Dissymmetry Theorem for Trees). Let T*{z) denote the expo¬ 
nential generating function of rooted trees (formed by taking the Euler derivative), 
T~{z) the collection of edge-rooted trees and T*~°{z) the exponential generating 
function of oriented edge rooted trees. We have the identity due to Otter m 

T'{z)+T-{z) = T{z)+T—°{z). (7.17) 

We realise that rooting a tree at an edge is the same as splitting the vertex 
set into two sets either side of the edge and having a rooted tree for each half i.e. 
T~ = £ 2 {T*) = ^'^ 2 ^ ■ realise that if we have an oriented edge then we have 
= {T'Y ■ This provides us with the modified identity: 

r%z)-\{T-{z)f=nz). (7.18) 

Corollary 7.7. We have the following relationship between doubly rooted and rooted 
trees: 

T*(z) 

T'%z) = -. (7.19) 

^ ’ 1-T*(z) ^ ’ 

Proof. If we take the Euler derivative of (I7.18p . then we obtain: 


r**(z) -T*(z)T** = r*(z), 

□ 

Theorem 7.8 (Functional Equation for Trees). Rooted tree generating functions 
satisfy a recursion relation 


T*( 0 ) = zexp(r*(z)). (7.21) 

In particular for the derivative of the tree generating function, we have: 

T'{z) = exp{zT'{z)). (7.22) 

This identity is found in [16]. 

Proof. Consider a rooted tree. The root may have an arbitrary number of neigh¬ 
bours. One can view the rest of the tree outside of the root as a collection of trees 
rooted at the corresponding neighbour of the initial root. Hence, we have the hrst 
z for the root of the tree and the exponential function since the order of the neigh¬ 
bours doesn’t matter so we need to divide by kl if we have k neighbours and for 
each neighbour we have a rooted tree. □ 

Corollary 7.9. We have that the second derivative of the tree generating function 
may be expressed in terms of the first: 

{T'{z)f 
1 - zT'{z) ■ 


T"{z) = 


(7.23) 
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Proof. If we differentiate ()7.22ll , then we obtain: 


T"{z) 

T"{z){l - zT\z)) 
T"{z) 


(T'(z) + zT"(z)) exp(zr'(z)) 
(T'(z))2 + zT'(z)T"(z), 

{T'iz)r 

l-zT'{z)' 


(7.24) 

□ 


Proof of Theorem [m We know that the sum of all tree graphs may be written in 
two ways using Cayley’s formula: 


n — 1 

^ ITpen.^WI =n"-2. (7.25) 

m—1 

We emphasise that any tree has a well dehned m for which is it m-splittable. 

We can express the left hand side easily in terms of coefficients of generating 
functions from (17.121) and (17.141) . For the right hand side we want the number of 
trees on n vertices to appear as the coefficient of This means we need to 

differentiate the tree generating function. Furthermore, the constant term of the 
left hand side is 0, but r'(0) = 1, so we must also take 1 and obtain: 

n—1 

(n - l)![z--i] ^ (Tpe„,i(z))- = (n - l)l[z-^]{r (z) - 1) (7.26) 

m—1 

which we have for every n. We also note that we may make the sum on the left 
hand side to infinity without affecting the outcome, since we will only have terms 
of greater powers of z from adding these factors. Hence we have the identity: 


2 pen,l(^) i 

1 ^Pen,l{z} 

This can be manipulated into the form: 

(7.27) 

Tpen.Az)-! 

If we differentiate (17.281). then we have the expression: 

(7.28) 

r (z) - 

We use corollary 17.91 to obtain: 

(7.29) 


(7.30) 

and using corollary 17.71 we have: 


TA^A^)=T-{z) + l. 

(7.31) 


Hence, we integrate both sides using the fact that Tpen,i(0) = 0 to obtain: 

Tpen,r(.) = ;^n"^^ 

(n + 1)! 

n—1 ^ ' 


(7.32) 
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□ 

7.3. Evaluating an Upper Bound. For stable potentials, we have the bound, 
from section 16.4[ 


l/3ra+l| 
(n + 1)! 


< 


m—1 



{TpenAz)r 


CAY 


[z”] 1 + -rpen,l{uz) 


(7.33) 


We may represent this coefficient through the contour integral representation: 

|."1 (l + irp..,,(„))” = J- £ ^ (l + iTp„,,(„«))” (7.34) 

where C is a contour around 0 within distance — from the origin, which is the radius 
of convergence of Tpen.iA^) ■ The right hand side can be bounded above by taking 
a contour with fixed radius R and since all coefficients are positive in Tpen,i(u2) 
the maximum value of the modulus along this circle is the function evaluated at R. 
We thus have: 


[Z"] 1 + -Tpen.l(uz) < 


1 + :^Tpen,l(w7?) 


R 


(7.35) 


for positive R < 

In order to find an optimal R, we differentiate and find the smallest positive root 
of the derivative equal to zero. That is we need to solve: 

(1 + iTpe„.i(ui?))—i(ui?) - 1 - iTpe„,l(ui?)) 


= 0 . 


(7.36) 


i?"+i 

Since iTpen,i(u7?) cannot equal —1, it follows that the equation is solved by: 

RTAr, i{uR) - 1 - iTpe„.i(ui7) = 0, 

’ u 

1 (ui?) - Tpen.l(ui?) = U. 

and hence we have the smallest positive solution to this equation dependent on u 
and call this c(u). This provides the upper bound for virial coefficients: 


lA 


'n+ll 


< 


C'(A)" /I + ^Tpen,l(uc(u)) 


(n+1)! (n + 1)! 


i(u) 


(7.37) 


This leads to the lower bound on the radius of convergence of the virial expansion: 

c(u) 


R-Wh > 




1 + irpen,l(uc(u)) 

8. Conclusions and Future Work 


(7.38) 


Combining the Penrose tree partition scheme with the algebraic relationship be¬ 
tween the virial and the cluster coefficients, we have been able to rigorously achieve 
Groeneveld’s lower bounds on the radius of convergence of the virial expansion for 
a classical gas with purely repulsive interactions. This may also be generalised to 
stable potentials through using a relabelling of the trees. This has also provided 
representation of the virial coefficients in terms of weighted trees. 
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Graph tree equalities, through the fundamental theorem of calculus of Brydges 
and Federbush [3], are an alternative way of understanding rewriting connected 
graphs in terms of trees. These have been generalised to matroids by Faris [5] and 
given a symmetric form by Abdesselam and Rivasseau [I]. These identities have 
not yet been made amenable to finding a product structure. A further possibility 
to deal with stable potentials is to use the ideas of Rivasseau [26] where a Kruskal 
algorithm gives an optimal way of writing a tree partition. 

Future work is to understand what one can achieve from other partition schemes 
and whether there is any sense of an optimal partition scheme from which one can 
obtain bounds. It is anticipated that future papers will understand the application 
of these to hardcore potentials in the discrete and continuous case, where the new 
tree weights should assist in making more accurate estimates. In addition, this new 
expression for virial coefficients in terms of trees should provide an easier method 
of making calculations than the method of Ree-Hoover bounds [5^ . 
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